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ABSTRACT 

Many intermediate age star clusters in the Magellanic Clouds present multiple main sequence 
turn-offs (MMSTO), which challenge the classical idea that star formation in such objects 
took place over short timescales. It has been recently suggested that the presence of fast ro- 
tators among main sequence stars could be the cause of such features dBastian & de Minkl 
120091) . hence relaxing the need for extended periods of star formation. In this letter, we com- 
pute evolutionary tracks and isochrones of models with and without rotation. We find that, 
for the same age and input physics, both kinds of models present turn-offs with an almost 
identical position in the colour-magnitude diagrams. As a consequence, a dispersion of ro- 
tational velocities in coeval ensembles of stars could not explain the presence of MMSTOs. 
We construct several synthetic colour-magnitude diagrams for the different kinds of tracks and 
combinations of them. The models that best reproduce the morphology of observed MMSTOs 
are clearly those assuming a significant spread in the stellar ages - as long as ~ 400 Myr - 
added to a moderate amount of convective core overshooting. Only these models produce the 
detailed "golf club" shape of observed MMSTOs. A spread in rotational velocities alone can- 
not do anything similar. We also discuss models involving a mixture of stars with and without 
overshooting, as an additional scenario to producing MMSTOs with coeval populations. We 
find that they produce turn-offs with a varying extension in the CMD direction perpendicular 
to the lower main sequence, which are clearly not present in observed MMSTOs. 
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1 INTRODUCTION 

There is now conclusive evidence that many star clusters in the 
Magellanic Clouds present double or multiple main sequence turn- 
offs (MMST O). The first hints o f thi s phenomeno n have been 
advanced by iBertelli et ail j2003h and iBaume et all d2007h . and 
were based on ground-based data for the LMC clusters NGC 2173 
and NGC 2154. Much more impressi ve and conclusive, how- 
ever, were the evidences brought by | Mackev & Brobv Nielsenl 
( l2007t) . lMackev et alj l20o3) . |Goudfrooii et all ( feOO^ for the clus- 
ters NGC 1846 NGC 1806, and NGC 1783 in the LMC, and 
iGlatt et alj j2008h for NGC 419 in the SMC, based on much deeper 
photometry obtained with the Advanced Camera for Surveys (ACS) 
onboard the Hubble Space Telescope (HST). In all these cases, the 
presence of a broad or MMS TO clearly stands o ut even from a vi- 
sual inspection of the CMDs. lMilone et al] j2009h revised the avail- 
able ACS photometry for 16 intermediate-age LMC clusters, find- 
ing evidence for broad turn-offs in about 70 % of them. 

The easiest interpretation of MMSTOs is that these clusters 
contain two or more generations of stars formed one after the other, 
over a time span of a few hundreds of Myr. This is indeed the inter- 



pretation adopted by the above-mentioned works. In support of this 
interpretation, there is the fact that a subsample of the intermediate - 
age clusters contain a dual clump of red giants l lGirardi et alj2009h . 
indicative of both the onset of electron degeneracy in stellar cores, 
and of a ~ 0.15 Mq spread in turn-off masses which is compati- 
ble with the one ind icated by the color spread in the turn-off region 
dRubele et alfeoiol. and work in preparation). 

This interpretation however is not an easy one, when one 
considers the challeng es it poses to the theory of star and clus- 
ter fo rmation (see e.g. iGoudfrooii et"al]|2009l ; IConrov & Spergell 
|2010|) . Note also that the prolonged star formation probably oc- 
curred in situ, within the relatively shallow potential well of < 
10 s Ma clusters, and are n ot due to the merging of different clus- 
ters dGoudfrooii et al]2009h . 

There is, however, an amazing coincidence regarding the age 
scales in all these cases. For all LMC clusters with the clear pres- 
ence of MMSTOs in the above-mentioned papers, the age interval 
between the youngest and oldest turn-off ranges from about 150 to 
~ 300 Myr, whe reas the clusters them selves have all ages between 
1 and 2 Gyr (see lMilone et~ai1l2009l) . NGC 419 in the SMC indi- 
cates a maximum age spread of 700 Myr, but again for a mean age 
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Figure 1. The evolutionary tracks computed for this work (dark solid lines), together with isochrones derived from them (narrow lines, in color in the 
electronic version), for the non-rotating, rotating, and overshooting cases (from left to right, respectively). The isochrones ages go, from top to bottom, from 
log(t/yr) = 8.5 to 9.3 at steps of 0.1 dex. As a reference to the eye, the isochrones with ages of 0.5, 1 and 2 Gyr are marked with a different colour. 



of 1 .5 Gyr dRubele et al.l201fj) . In terms of CMD features, the MM- 
STOs always cover a magnitude range of about ~ 0.5 mag in pass- 
bands like V or /, and ~ 0.15 mag in colours like V—I. The MM- 
STO phenomenum is observed in star clusters located overall over 
the LMC disk and also in the SMC, although i t seems to be pref- 
erenti ally observed in the most massive objects dConrov & Spergell 
The homogeneity in the mean ages of clusters with MM- 
STOs could be indicating dramatic large-scale events in the past 
history of the Magellanic Clouds, that could have affected the for- 
mation of many of its more massive clusters in a similar way. Or, 
perhaps, this coincidence is simply suggesting us to look for a com- 
plet ely different interpre t ation t o the MMSTO phenomenum. 

iBastian & de Minkl 1 12003) have so far advanced the only al- 
ternative explanation for such CMD features, namely that the MM- 
STOs appear in aggregates of coeval stars, due to the color disper- 
sion caused by the presence of fast rotators among the main se- 
quence stars. Their suggestion is based on simple considerations 
about the shape of stellar evolutionary tracks with and without ro- 
tation, on some observational facts, and on some less well-justified 
assumptions, namely that: 

• a significant fraction of main sequence stars rotate with veloc- 
ities of about ~ .4 times the critical break up velocity, fi cr it ( see 
iRover et al.l2007l) ; 

• rotation causes stars to shift, in general, to redder colours and 
slightly brighter luminosities in CMDs; 

• for the most massive stars, close to the MSTO, the effect of 
rotation is quite similar throughout the main-sequence evolution, 
so that the effects are considered as displayed when stars reach the 
end of the main sequence; 

• the effect of rotation gets less important at smaller masses, and 
disappears for masses smaller than ~ 1.2 A4©. 

Note that the combination of all these different effects is required to 
explain the shape of the observed MMSTOs. In addition, strongly 
bi-modal distributions of rotational velocities are needed to explain 
the few cases in which the turn-off appear as a double feature. 

In this paper, we investi gate whether rotation doe s really be- 
have in the way sugg ested bv lBastian & de Minkl d2009h . Our anal- 
ysis is based on isochrones computed with detailed models of ro- 
tating stars (Sect. |2), and do not suffer from the approximations 
adopted by them. In Sect. [3] we compare the performance of co- 



eval models with a dispersion of rotation, with that of non rotating 
models with a dispersion of ages, in producing the MMSTO feature 
observed in the LMC cluster NGC 1846. The basic conclusions are 
discussed in Sect. [4] 



2 MODELS OF ROTATING STARS AND ISOCHRONES 

We start computing a series of evolutionary tracks with initial 
masses in the interval relevant for MMSTOs. An initial heavy ele- 
ments mass fraction of Z = 0.008, suitable for LMC populations 
of young-to-intermediate ages, is adopted. 

The main-sequence evolution of models with and without ro- 
tation is computed by u sing the Geneva stellar evolution code 
feggenberge r et al. 2008). The stan dard mixing-length formalism 
for convection jBohm-Vitensei 19581) is used with a solar calibrated 
value of the mixing-length parameter. Overshooting from the con- 
vective core into the surrounding radiatively stable layers is not 
included. Rotating models are computed with an initial velocity on 
the ZAMS of 150 km s - , which corresponds to a typical value for 
rotating s tars in the mass in terval between 1.5 and 3 Mq studied 
here (e.g. IRover et al.ll2007h - Fl Note that the ratio u = Q,/Q, CI \ t 
of these models remains always inferior to 0.7, so that the effects 
related to the variation with colatitu de of the effective temperature 
are negligible tekstrom et ai1l2008h . These models include a com- 
prehensive treatment of shellular rotation (Zahn 1992) and merid- 
ional circulation is treated as a truly adv ective process (for more 
details, see e.g. lEggenberger et al.ll2O10h . The rotating and non- 
rotating models share the same initial parameters except for the 
inclusion of shellular rotation. 

These two sets are complemented by one without rotation but 
with overshooting from the convective core, computed with the 
same code. The overs hooting parameter is set t o a ov = 0.25 pres- 
sure scale heights (see lMaeder & Mevnetll 19891 for details). 

The tracks are presented in Fig.[T] They cover all evolutionary 



1 Note however that for the lower metallicities of the LMC and SMC, the 
median rotation velocit y might be somewh at higher than in the Solar Neigh- 
bourhood sample from lRover et all fe007l) . This difference however is un- 
likely to affect the main arguments discussed in the following. 
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Figure 2. Left panel: HR diagram of the different isochrones obtai ned in Fig.[T] for the same age of 1.58 Gyr. Right: a synthetic HR diagram where the same 
isochrones are populated by 10 4 M Q following a lChabrieil feOOll) log-normal IMF. To avoid the excessive superposition of points in the diagram, they were 
displaced by random Gaussian errors of 0.01 dex in log L and 0.0025 dex in log T c ff . 



stages from the zero-age main sequence (ZAMS) up to either He- 
ignition, or the RGB-bump of less massive stars. 

The same figure presents a set of isochrones derived from 
them. Despite the relatively large mass interval bet ween succes- 
sive t racks, the interpolation algorithm (the same as in lGirardi et al.l 
l200d) works very well and produces isochrones with all the fine 
details of the original tracks, especially along the MSTO region. In 
order to include the lower main sequence in th e isochrones, the new 
grids of tracks have been complemented with lGirardi et alj J2000h 
non-rotating tracks of M ^ 1.0 Mq, which, for the ages consid- 
ered in this work, are always very close to the ZAMS. 

The two sets of tracks without overshooting make the basic 
set with which we will discuss the effect of rotation on the CMDs. 
We recall that these tracks do follow, qualitatively, the changes 
in the stellar T e s and L due to rotation that are at the base of 
iBastian & de Minkl J2009h work. 



3 SYNTHETIC CMDS AND COMPARISON WITH DATA 

Let us now consider the effect that rotation has in the CMDs of star 
clusters. Let us assume for the moment that all the stars in clusters 
with multiple turn-offs have exactly the same age and initial chem- 
ical composition. The theoretical entities suitable to describe such 
a situation are the stellar isochrones. In the following, we will cen- 
ter our discussion on isochrones computed for an age of 1.58 Gyr, 
however we advance that very similar conclusions would be ob- 
tained for any age between 1 and 2 Gyr. 

In the left panel of Fig. [2] we superpose isochrones of the 
same age computed for the 3 cases mentioned in Sect. [2] in the 
HR diagram. The result is somewhat su rprising: contrarily to what 
sustained in IBastian & de Minkl ( |2009|) . models with rotation do 
not present a cooler and fainter turn-of^ with respect to their 
non-rotating counterparts. Instead, they have a slightly hotter and 
brighter turn-off. This result cannot be easily inferred while look- 
ing only at the evolutionary tracks in Fig. Q] which at first sight 
seem to suggest that rotating models have a broader main sequence, 
hence being able to exten d to cooler temperatures, as sustained in 
IBastian & de Minkl i2009h . This is true only for pairs of evolution- 
ary tracks of the same mass. What we see in the isochrones is dif- 



2 Throughout this paper, turn-off is intended as the tip of the locus contin- 
uously populated by stars along the main sequence. 



ferent, it is the combined effect of changes in the morphology of the 
evolutionary tracks, and of the changes in their lifetimes. 

The increase in the main sequence life times caused by ro- 
tation, in particular, is q uite remarkable (see IPalacios et al]|2003l : 
lEggenberger et al]|2010h and constitutes one of the main factors in 
play here. As discussed in lEggenberger et al] ( l2010h . classical mod- 
els with a v ro t — 150 km/s have a similar increase in their main se- 
quence lifetimes as models with a moderate amount of overshoot- 
ing, that is a ov = 0.1. Su ch an increase in the lifetim es of rotating 
stars has been ignored in IBastian & de Minkl J2009h . and this has 
been determinant in leading to their different conclusions. 

In the right panel of Fig. [2] we present synthetic HR diagrams 
corresponding to 1.58-Gyr old clusters with and without rotation 
- in this latter case assuming that all stars rotate with the same 
initial velocit y - and wi t h over shooting, for a total mass of 10 4 Mq 
following the lChabrierl j200lh log-normal IMF. This panel clearly 
evinces what are the regions of the diagram effectively populated 
by stars. In particular, one may notice that the differences between 
non-rotating and rotating models is not dramatic, if one considers 
only the region where the MSTO is drawn. 

Then, using the TRILEGAL code dGirardi et ai] 120051) . we 
simulate populations in the CMD in which MMSTOs are usually 
observed, namely the F555W vs. F555W — F814W (approximately 
V vs. V — /) diagram for the filters in the the ACS/WFC camera on- 
board HST. Such simulations are shown in Fig. [3] The bolom etric 
correction tables in use are described in lGirardi et ai] d2008h . The 
simulated photometry is displaced by the absolute distance modu- 
lus of 18.5 mag and a l/-band foreground extinction 0.2 mag, which 
are values typical of LMC star clusters. Then, photometric errors of 
o = 0.01 mag are applied to the photometry, so as to help in the 
visualisation of the densest parts of the CMDs without however 
modifying their morphology in any appreciable way. All simula- 
tions include a fraction of 20 % of binaries with mass ratios uni- 
formly distributed between 0.4 a nd 1. The t otal m ass is assumed to 
be of 2.4 x 10 4 M Q , following a lChabried J200lh log-normal IMF. 

The models in the top row of Fig. [3] from (a) to (c), are 
for a single age of 1.58 Gyr, and using each one of the 3 sets 
of tracks presented in this work. These top panels illustrate the 
modest effect that binaries have in these CMDs, in a way similar 
to the synthetic CMDs already presented by many authou rs (e.g. 
iBertelli et alj2003l ; iKerber et alj|2007l ; iMackev et alj2008l) . 

Panels in the second row of Fig. [3] represent attempts to pro- 
duce MMSTOs by adding stellar models with different properties 
but again for the same single age. Panel (d) shows the case of mod- 
els with and without rotation and otherwise identical physics. This 
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Figure 3. Panels (a)-(i) present synthetic CMDs illustrating the effect of different prescriptions regarding the distributions of initial rotational velocities v TO t , 
overshooting efficiency a ov , and ages t among the clusters stars. All these models assume a modest fraction of binaries and small p hotometric errors. For 
comparison, panel (j) shows the ACS/WFC data for the central region of the LMC star cluster NGC 1846 (from lGoudfrooii et alj2009l) , with its well defined 
MMSTO. For each panel, we provide the number of objects in the turn-off region defined by 19 < F814W < 21.5 and 0.2 < F555 - F814W < 0.8. The 
top row shows models for a single age and for (a) no rotation without overshooting, (b) rotation without overshooting, and (c) no rotation without overshooting. 
The second row show mixed models, again for a single age: a mixture of (d) rotating and non rotating stars, without overshooting, (e) stars with and without 
overshooting, without rotation, and (f) non-rotating stars with overshooting, together with rotating stars without overshooting. It is evident that none of these 
cases produces the peculiar "golf club" shape of the observed MMSTO. In the third row, from (g) to (i), we have the same models as un the top row, but now 
with the star formation spanning and interval of 400 Myr. The similarity with the observed MMSTO is evident, especially in the panel (i), corresponding to 
models with a moderate efficiency of convective overshooting but no rotation. 
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corresponds to the situation favoured bv lBastian & de Minkl fioogj) 
as being at the origin of MMSTOs. As can be noticed, our results 
produce a dramatically narrow MSTO, clearly very different from 
both t he observed one and those illustrated in lBastian & de Minkl 
d2009h . 

Panel (e) instead explores the case in which there is a mix- 
ture of stars with different efficiencies of core overshooting, but 
observed at the same age. In this case the MSTO opens into a 
composite feature, with two distinct MSTOs. However, the brighter 
MSTO (corresponding to the overshooting models) is clearly much 
more extended to larger radii (towards the top-right corner of the 
figure) than the fainter MSTO. This again contrasts with observa- 
tions, in which the MMSTOs are seen to have about the same ex- 
tension in the direction perpendicular to the lower main sequence. 
Finally, the panel (f) shows the combination of non-rotating stars 
with overshooting, together with rotating stars without overshoot- 
ing, so as to simulate an hypothetical reduction of the extent of 
the overshooting region driven by rotation. This combination pro- 
duces CMDs slightly more similar to the observed one, but does 
not reproduce its detailed shape; in particular, it can be noticed that 
the fainter turn-off, corresponding to rotating models without over- 
shooting, is again less extended in the direction perpendicular to 
the main sequence than the brighter one. This discrepancy cannot 
be remediated imposing that the fast rotators have some residual 
amount of overshooting, since this would increase their lifetimes 
and move the fainter turn-off towards the brighter one, creating a 
situation similar to the one in panel (d). The adoption of faster rota- 
tional velocities in models without overshooting would also move 
the faintest turn-off upwards in a similar way. 

The third row of Fig. [3] from (g) to (i), presents the models 
derived from each set of homogeneous tracks, but now assuming 
that star formation proceeded at a constant rate in the age interval 
from 1.3 to 1.7 Myr. The similarity with the observed MMSTO is 
now evident, especially in the panel (i), corresponding to models 
with a single efficiency of overshooting but no rotation. 



4 FINAL CONSIDERATIONS 

The present work demonstrates that stellar rotation, as predicted by 
detailed evolutionary models, produces featur es in CMDs which 
are qu ite different from those illustrated in iBastian & de Minkl 
J2009h . We find that, despite the changes that rotation causes in 
the HR diagram of evolutionary tracks, their longer lifetimes make 
the isochrones derived from these tracks to have a MSTO posi- 
tion almost indistinguishable from those derived from non-rotating 
models. It follows that coeval isochrones derived from rotating plus 
non-rotating tracks - with otherwise identical physics - are not able 
to describe the detailed shape of the observed MMSTOs in Magel- 
lanic Cloud clusters. 

On the other hand, Fig. [3]an d the careful CMD-r econstruction 
work bv lRubele etail ( 1201 Ol) and lRubele etail J201 J) demonstrate 
that the other possible explanation for the origin of MMSTOs - 
namely prolonged histories of star formation taking place in the 
clusters - does provide excellent qualitative and quantitative fits of 
the CMD in the clusters NGC 419 and NGC 1751. This despite 
these latter clusters present dual red clumps in addition to MM- 
STOs, which represents a big additional challenge to CMD-fitting 
algorithms. It is hardly conceivable, after the results from this work, 
that a coeval population of stars with a dispersion in their rotational 
velocities, could produce a CMD-fitting that even approaches the 
quality reached in those CMD-fitting works. We are left with no 



better alternative, for the moment, than accepting that prolonged 
star formation has probably occurred in these clusters, although 
we are aware of the great challenges it represents for dynamical 
models of stellar and cluster formation. 

We are also aware that stellar rotation - and fast rotators - 
occurs anyway among the cluster stars. What is evinced by our 
models is simply that the effect of fast rotation on CMDs is mod- 
est and more likely similar to the colour and magnitude dispersion 
caused by bi naries, rather than the dram atic drawing of MMSTOs 
advocated bv lBastian & de Minkl J20091) . Following our reasoning, 
the detection of fa st rotators at the MMSTO s is just expected, and 
would not suvvort \Bastian & de M ink 2009 ' hypothesis unless the 
fast rotators are all observed at the reddest MSTO. This is an as- 
pect to be tackled by future obse rvations of rotation i n cluster stars, 
as those recently announced bv lPlatais etail j201Ch for the open 
cluster Tr 20. 

In conclusion, the explanation based on prolonged periods of 
SFH seems to be the only surviving one, in this moment, to ex- 
plain the origin of MMSTOs in Magellanic Cloud clusters. De- 
spite the perplexities and difficulties it may cause to star for- 
mati on theories, the pre sence of dual red clumps in some clus- 
ters l lGirardi et alJl2009T) and the high qualit y of the resu l ts ob- 
tained from the m odelling of their CMDs dRubele et all l20ld : 
iRubele et alj201 lb strongly favour this interpretation. Detection of 
chemical self-enrichment in clusters with MMSTOs would likely 
be the decisive data to confirm or disprove this explanation. 
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